Electron penetration in the nucleus and its effect on the quadrupole interaction 
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A series expansion of the interaction between a nucleus and its surrounding electron distribution 
provides terms that are well-known in the study of hyperfine interactions: the familiar quadrupole 
interaction and the less familiar hexadecapole interaction. If the penetration of electrons into the 
nucleus is taken into account, various corrections to these multipole interactions appear. The best 
known one is a scalar correction related to the isotope shift and the isomer shift. This paper discusses 
a related tensor correction, which modifies the quadrupole interaction if electrons penetrate the 
nucleus: the quadrupole shift. We describe the mathematical formalism and provide first-principles 
calculations of the quadrupole shift for a large set of solids. Fully relativistic calculations that 
explicitly take a finite nucleus into account turn out to be mandatory. Our analysis shows that 
the quadrupole shift becomes appreciably large for heavy elements. Implications for experimental 
high-precision studies of quadrupole interactions and quadrupole moment ratios are discussed. A 
literature review of other small quadrupole-like effects is presented as well. 

PACS numbers: 



I. INTRODUCTION 

Atomic nuclei are no mathematical point charges, but 
objects with a shape and a size. This affects the way 
in which they interact with electrons, especially when 
electrons penetrate the nuclear volume and render the 
usual 'far-field' approximation invalid. These 'near-field' 
effects lead to tiny corrections to all terms in the multi- 
pole expansion for the electrostatic interaction between 
nuclei and electrons. The correction to the monopolc 
term corresponds to experimentally well-known phenom- 
ena: the isotope shift in atomic spectroscopy and the 
isomer shift in Mossbauer spectroscopy. An analogous 
correction to the quadrupole term - coined here the 
quadrupole shift [96[ (QS) - should exist as well. The 
existence of such an effect has been touched upon a few 
times in the literature of the past decades p], 0, H, 0, H, H] ; 
but to our knowledge a systematic study is lacking. In 
this paper, we present a mathematical treatment of the 
quadrupole shift by a twofold application of first-order 
perturbation theory, which leads to a simple analytical 
expression (Sees. HHlIII[) . We point out that in order to 
compute numerical values for the quadrupole shift from 
first-principles, it is necessary to perform fully relativis- 
tic calculations that take explicitly a finite nucleus into 
account (Sec. IIV|) . Density Functional Theory calcula- 
tions of the quadrupole shift for a set of simple crys- 
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tals show that the size of the quadrupole shift strongly 
grows with the mass of the isotope, an effect that turns 
out to have an electronic rather than a nuclear origin 
(Sec. |V}. Except for the heaviest elements (actinides), 
the quadrupole shift is only a minor correction to the 
quadrupole interaction. We discuss how it shows up 
in experiments, and how it could possibly be exploited 
to improve the accuracy of experimentally determined 
quadrupole moments (Scc. lVl|) . Especially for the exper- 
imental determination of ratios between nuclear electric 
quadrupole moments, the accuracy that can be reached 
by the most precise molecular beam spectroscopy ex- 
periments is good enough to make it relevant in some 
cases to take quadrupole shift corrections into account. 
The quadrupole shift is only one of a set of small effects 
that can affect the regular quadrupole interaction. The 
(sometimes fairly old) literature on these other effects is 
reviewed in App. [A] We suggest that for high-precision 
studies it is relevant to revisit these small quadrupole-like 
perturbations with modern computational methods. 



II. FORMALISM 

A. Classical interaction energy without 
charge-charge overlap 

The classical electrostatic interaction energy between 
a positive (nuclear) charge distribution p(f) and a poten- 
tial v(r) due to a surrounding (electron) charge distribu- 
tion n(F') is formally given by (with eo being the electric 
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constant) 

E = I p(f)v(f)dr 



47re 



«g*£W, (i) 



and can be expressed by the standard multipole expan- 
sion in spherical harmonics Q: 



(2) 



with r< = min(r, r') and r> = max(r, r'). This leads to 
an infinite sum of double integrals, each with the dimen- 



sion of energy: 



E — E-2i — Eq + Ei + E4 
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(Odd terms will vanish in the cases of interest here, see 
Sec. Ill Bl ) It is the second term E2 that will be of interest 
in the present work: 



Eo = hvn = 
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The frequency vq is experimentally accessible, and is 
called the nuclear quadrupole coupling constant (NQCC). 
Due to the varying assignment of r< and r> to 'nuclear' 
(r) or 'electron' (r') coordinates, quantities as E2 are an 
intricate mixture of properties of both charge distribu- 
tions p(r) and n(f). Only in the special case where both 
charge distributions do not overlap (r< = r and r> = r '), 
Eq. ([I]) can be written in terms of properties that depend 
entirely on only one of the charge distributions: 
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where Qi m and Vi m are the components of the nuclear 
multipole moment and electric multipole field tensors of 
rank I, respectively: 
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r l p(r)Y lm (n)dr 
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T n(f')^m(fi')<8"- (7) 



When this formalism is applied to describe nuclei and 
electrons, the simplification by Eq. ([5]) can never be 

made: s-electrons and relativistic pi -electrons have a 

2 

non-zero probability to appear at r = 0, and therefore the 
nuclear and electron charge distributions always overlap. 
Nevertheless, motivated by the very small size of the re- 
gion where this overlap happens compared to the volume 
of the rest of the atom, one can in a first approxima- 
tion neglect this concern and apply Eq. ([5]) to atoms, 
molecules and solids. This is where the concept orig- 
inates of an electric- field gradient (EFG) tensor (Vz m ) 
that interacts with a nuclear quadrupole moment tensor 
(Qsm) to produce an experimentally observable interac- 
tion energy (E2). Although E2 itself is a well-defined 



p(f)n(f')^Y* (Q)Y 2 , m (n')dfdr'. 



(4) 
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observable property, its description by a quadrupole in- 
teraction energy only 



+2 
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rather than by Eq. ((J) is an approximation. 



(8) 



B. Overlap corrections 

We will now derive explicit expressions for the correc- 
tions that need to be added to Eq. ([5]) to obtain Eq. ^ 
(and similarly for other values of I). Rather than us- 
ing the multipole expansion in spherical harmonics from 
Eq. ([2]), we start from a Taylor expansion of the electro- 
static potential v(r) = l/(47reo) / n(f')/\r — r'\ dr' in 
the interaction energy of Eq. JT]): 



E = J p(r)v(r)dr = v(0) J p(r)df 
+ ^Uj(0) / Xip{f)dr 

■ ^rX!'V :l)i / x l x J p(f)dr 
+77y y^%fc(Q) J XiXjX k p(r)dr 

i,j,k 

+i ^2 J XiXjX k xip(f)dr + 0(6). (9) 

i,j,k,l 

In order to recognize in this expression the multipole mo- 
ments and multipole fields from Eq. ^ and Eq. (JT]) , one 
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has to make substitutions like this one (the example is 
for the quadrupole moment): 

1 



XiXjp(r)dr 



{ZxiXj - r 8 lj )p{r)dr 



+ o / r 2 p(r)dr 8 i:j , 



(10) 



where Qij are the components of the quadrupole tensor 
Qim (Eq. ([6])), but now in Cartesian form. This yields 
for the first three even orders in Eq. © the following 
nuclear multipole moments in Cartesian form: 



M 



H 



ijkl 



p(r)v(r)df = eZ (11) 

(ZxiXi - r 2 5 l j)p{r)dr (12) 
3 • 5(7 Xi 

-f H {xi,Xj,Xk,Xi))p{r)dr (13) 



with f H (xi,Xj,x k ,xi) = r 2 

XjXkSti + XjXiSik + XkXiSij 
The corresponding electric multipole fields in Cartesian 
form are: 



XiXjSki + XiXkSji + XiXiS k: j + 



SijSki + SikSji + SuSjk 



V = v(0) (14) 
Vij = (didjv{0) -^ASij) Av(0) (15) 

Vijkl = didjdkd^-f^AviO) (16) 



with fYjM = 

djdAk+dkdAj 



didjSki + d l d k S J i + didiS kj + djd k Sii + 

The exprcs- 



SijSki+SikSji+SuSjk 
sions in Eqs. (IT]) to jl3]) and Eqs. (14]) to (16]) arc iden- 
tical to the ones in Eqs. ([6]) and (7]), respectively. They 
have the same number of degrees of freedom: 1, 5 and 9 
for the zeroth, second and fourth order moment/field. 

After having inserted into Eq. ([9]) all substitutions as 
in Eq. (|10p . the interaction energy can be written as 

E = M^V + ^{r 2 }Av{0) + ^{r 4 }A 2 v(0) 
MI 



MS^ 1 
^3 QijVij 



MS< 2 



Qi 



i ^{(xiXj - ^5 l3 )r 2 md 3 - f %)A«(0) 
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0(6), 



(17) 



HDI 



where all integrations over the nuclear charge density p(r) 
are noted in short-hand by {curled brackets}. Eq. (TTT|) 
contains no odd order terms (dipole, octupole,. . .), since 
nuclei have no odd order electric moments due to time 
reversal symmetry 0]. We see that Eq. p7|) contains 
dot products between multipole moments and fields as 
in Eq. ([5]): the monopole (MI), quadrupole (QI), hexade- 
capolc (HDI),. . . interactions. These are the only contri- 
butions in the case without charge-charge overlap. Ad- 
ditionally, an infinite set of even order correction terms 
appears now as well - due to parity, there are no odd 
order corrections. In Tab. [H a general naming system 
and a corresponding set of symbols are presented: the 
n th order quasi multipole moment multiplied (dot prod- 
uct) with the n th order quasi multipole field leads to the 
n th order multipole shift. From the general trends in this 
table one can infer the structure of the higher order cor- 
rections that were not explicitly derived in Eq. (|17[) 
they are shown in the table in red. 

There is a qualitative difference between the multipole 
fields in the first column of Tab. [I] and the quasi mul- 
tipole fields in all other columns. The multipole fields 
depend on the potential i;(0) at the nucleus, which de- 
pends via integration on the charge distribution every- 
where else in the system. Multipole fields are therefore 
integrated quantities, determined by the entire density. 
The quasi multipole fields depend on the Laplacian of the 
potential at the nucleus (Aw(0)), which is by the Poisson 
equation proportional to the electron charge density at 
the nucleus (n(0)) (Au(0) = — n(0)/eo). Quasi multipole 
fields are therefore point quantities, determined by the 
electron density in a single point only. 

In the next section, the results of Eq. (fT?)) and Tab. U 
for a system of two classical charge distributions will be 
translated to a quantum formulation. This will make it 
applicable to atoms, molecules and solids. Known exper- 
imental consequences of these overlap correction terms 
will be summarized in Sec. IIIIl The core of the present 
work deals with the first order quadrupole shift QS^\ 
which is the first order correction to the quadrupole in- 
teraction. 



C. Quantum formulation 

In order to translate Eq. (fTTJ) to quantum mechan- 
ics, Hamiltonian operators corresponding to all its terms 
are required. The structure of Eq. (fT?]) suggests a 
perturbation theory treatment, with the monopole in- 
teraction term as the unperturbed Hamiltonian, and 
the other terms as small perturbations. The monopole 
term depends via r° on the (small) nuclear coordinate 
(r cx 10~ 15 m) and via 1/r' on the electronic coordinate 
(r' oc 10~ 10 m). Among all small corrections in Tab. [H 
the two largest ones are the quadrupole interaction QI 
and the first order monopole shift MS' 1 - 1 - both have a r 2 
in their nuclear parts and a second derivative of the elec- 
trostatic potential (leading to 1/r' 3 ) in their electronic 
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TABLE I: Systematic overview of nuclear multipole and quasi multipole moments and electric multipole and quasi multipole 
fields that appear in the multipole expansion of two interacting (and overlapping) classical charge distributions. The first 
column gives the regular multipole expansion for point nuclei: the monopole, quadrupole and hexadecapole interactions. The 
next columns give the quasi multipole moments/fields for every multipole interaction, denoted by a tilde: these are corrections 
to the multipole interactions due to electron penetration into an extended nucleus. Colored text is by generalization only, and 
is not systematically derived in this work. The objects in each line are spherical tensors of a given rank (rank for line 1, rank 
2 for line 2, rank 4 for line 3, . . .). 



Order 


Multipole 
moment 
/ field 


First order 
quasi moment 
/ quasi field 


Second order 
quasi moment 
/ quasi field 


O(0) 


MI: 

M oc r°y 00 
V ocu(O) 


MS (1) : 
M« oc {r 2 Y 00 } 
V" (1) oc Au(0) 


MS (2) : 
M< a > oc {r 4 Y 00 } 
t/ (2) oc A 2 v(0) 


0(2) 


QI: 

Q oc r 2 Y 20 
Vij oc dijv(0) 


QS«: 
Q (1) oc {r 4 Y 2Q } 
V<-p oc d tJ Av(0) 


Q {2) oc {r 6 y 20 } 
V} 2) oc dij A 2 v(0) 


0(4) 


HDI: 

H oc r 4 Y 40 
Vijki oc d ijk iv(0) 


HDS (1) : 
H w oc {r 6 Y 40 } 
« d ijhl Av(0) 


HDS (2) : 
H {2) oc {r 8 Y 40 } 
oc d ljk ,A 2 v(0) 



parts. These two leading corrections will be taken as the 
small perturbation. 

The Hamiltonians that correspond to the entries in 
Tab. |T] operate on the direct product space of wave func- 
tions for the nuclear and the electron subspaces. The 
ground state of the monopole Hamiltonian is a direct 
product between the nuclear ground state and the elec- 
tronic ground state wave function. With M = eZi 
(Eq. pip . 11 is the identity operator on the nuclear space) 
and V = v(0) (Eq. (fT4")) . v(0) is an operator on the elec- 
tronic space that returns the potential at f= due to a 
given wave function 'J'), the unperturbed monopole in- 
teraction Hamiltonian is 

Hmi = eZ 1®«(0). (18) 

Evaluating this for the ground state wave function 
\I (g> ^o) of the combined nuclear+electronic system leads 
to (|7) is the ground state of the nucleus, and |^o) tnc 
ground state of the electron system with a point nucleus): 

= (I\eZ n|i)-<*o|fi(0)|*o> 

= eZv(0) , (19) 

which is the leading term in Eqs. or (fT7|) . The label pn 
('point nucleus') emphasizes the difference with Eq from 



Eq. ([3]). The quantity v(0) - the electrostatic potential 
at the nuclear site for a point nucleus - is accessible by 
first-principles codes. 

The perturbing Hamiltonian is (see Tab. [J for the no- 
tation) : 

Hp = Hqi +W MS (i). (20) 

In first order perturbation theory, the energy corrections 
due to this perturbation are found by evaluating the per- 
turbing Hamiltonian in the ground state of the unper- 
turbed Hamiltonian. Assuming a non-degenerate ground 
state in the electron subspace, it is advantageous to write 
the Hamiltonians immediately in a more familiar form 
where the electronic matrix elements are already evalu- 
ated and are treated as known (^computable) quantities. 
After similar algebra as for the monopole Hamiltonian, 
this leads to this form for the monopole shift Hamiltonian 
(it contains the mean square radius (r 2 ) of the nucleus 
and the electron density n(0) at the position of the nu- 
cleus): 

n Msm = ~n(0) (r 2 )i (21) 

The quadrupole Hamiltonian TLqi contains the (spectro- 
scopic) quadrupole moment of the nucleus Q and the 
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quadrupole field of the electrons V zz (principle compo- 
nent of the electric-field gradient tensor) (see e.g. Ref.Q): 



eQV Zi 



4(2/ - l)Ih z 



(22) 



Diagonalizing these two Hamiltonians in the nuclear 
states leads to the desired energy corrections in first order 
perturbation. Formally: 



= E'P n + (I\H 



MS* 1 ) 



■Hqi\I) 



= £T + {I\H Msm \I) + {I\H QI \I) 



E 



E, 



[i] 
Qi- 



(23) 



Here, E^j s[1) is a correction to the monopolc energy E^ n 
for a point nucleus due to (s- or pi-)electron penetration 
into the volume of a spherical nucleus. The quadrupole 
interaction energy Eq\ is a correction due to the devia- 
tion from spherical symmetry of this nucleus. 

There is a second group of entries with even much 
smaller corrections in Tab. [Q the HDI, QS^ and MS^ 2 ) 
terms all have r 4 and 4 derivatives of the electrostatic 
potential (— ► 1/r' 5 ). The corresponding Hamiltonians 



7~Lhdi 



128/(7- 1)(2/- 1)(27 — 3)?l 
1 eQn,, 



14e 4(2/ - l)ITf 
-||An(0) (r 4 )ll. 



35/ 4 - mill 2 + 3/ 4 + 25h 2 I 2 - 6h 2 I 2 



(si 2 i 2 ) 



(24) 

(25) 
(26) 



r 



The (diagonal part of the) hexadecapole Hamiltonian, 
Eq. ([24]) . is taken from the literature the quadrupole 
shift Hamiltonian, Eq. (|2"5")l . is derived explicitly in [i| 
and similar algebra as for the first order monopolc shift 
Hamiltonian leads to the second order monopolc shift 
Hamiltonian, Eq. (|26[) . As they are much smaller than 
the QI and MS^ 1 - 1 terms, it makes little sense to add 
these corrections to the Hamiltonian of Eq. (f2T)|) right 
away. Rather one should consider a first order pertur- 
bation to the Hamiltonian of Eq. (T2"D|) , which itself was 
already a perturbation to the monopole Hamiltonian of 
Eq. (JTHJ) . This means: find the perturbed eigen states of 
Eq. (|20p in first order, and evaluate the new perturba- 
tions as given by the Hamiltonians in Eqs. P4"|) - (f2"6")) in 
these eigen states. In the present work, we are interested 
in the first place in H.q S {i) , as it has the symmetry of 
a quadrupole interaction: this Hamiltonian, evaluated in 
the ( approximate ) eigenstates for a system with a finite 
and quadrupolarly deformed nucleus, gives an additional 
contribution to the regular quadrupole interaction. It can 
be interpreted as the influence on the quadrupole inter- 
action of electron penetration into the nuclear volume: 
the quadrupole shift. The quadrupole shift Hamiltonian 
of Eq. (|2"5|) expresses the influence of the finite nucleus 
on the multipole expansion. Evaluating this Hamiltonian 
for a density obtained from a first-principles calculation 
with a finite nucleus expresses the influence of the finite 
nucleus on the electronic wave functions. 

There is an alternative way to express this same effect: 
consider the Hamiltonian of Eq. (|2"D|) up to second order 
perturbation. Among others, the second order energy ex- 



pression will contain a cross term between QI and MS' 1 ), 
which has the same symmetry as the quadrupole interac- 
tion (this can be easily seen because the Yoo term of the 
monopolc shift is a scalar quantity that does not change 
the symmetry). Compared to the previous strategy this 
method has the advantage that the same Hamiltonian 
is kept, but the disadvantage that second order matrix 
elements in excited states have to be evaluated. It is 
technically easier to evaluate a new perturbation in the 
ground state of the previous perturbation. The underly- 
ing physics, however, is the same. 



The second order perturbation description has been 
applied in 1970 by P. Pyykko for approximate and non- 
relativistic calculations in a few test molecules (see also 
Fig. HI) [l|. The first order + first order perturbation 
description has been used in 2003 by Thyssen et al. || 
for the case of Lil, albeit in an implicit way that did 
not clearly showed the twofold application of first order 
perturbation theory The twofold application of pertur- 
bation theory will be the method used in the present 
work as well, not at least because it leads to concise ana- 
lytical formulas. In 2006, also Karl and Novikov derived 
the so-called "contact terms" of the quadrupole interac- 
tion. They used the Feynman diagram technique and 
evaluated the results for hyperonic atoms [l(| El ■ Our 
derivation was made completely independent from the 
ones by Thyssen et al. and Karl and Novikov, and the 
observation that the final expressions agree is a strong 
test of mutual correctness. 
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D. Zooming in on E2 

The regular quadrupole interaction and first order 
quadrupole shift together provide our approximation to 
E 2 : 

E 2 = hv Q w Eqi + E QS u) = hvQi + hv QS m. 

(27) 

Both terms consist each of a product between a nuclear 
quantity and an electron quantity. As this shows to which 
nuclear and/or electronic properties one gets access by 
measuring E2, wc discuss them now. The two relevant 
nuclear quantities are (see Eqs. (|2"2"|) and l[2"5"j)): 



'mi = ±; 



Hqi -> eQ = I p(r)(3z 2 - r 2 )dr oc (r 2 Y 20 ) 

(28) 

p(r)(3z 2 - r 2 )r 2 dr oc (r 4 Y 20 ). 

(29) 



The quasi quadrupole moment Q has an additional r 2 
in the integral compared to the quadrupole moment 
Q. It is therefore a quantity that bears similarity 
with the quadrupole moment (r Y20) (through the In- 
dependence) as well as with the hexadecapole moment 
(r 4 Y4o) (through the r 4 -dependence). 

The corresponding electronic quantities are: 



Kqi 



(30) 



= ~ (d zz - j ) 



The integrated quantity (cf. Sec. Ill B|) V zz is the princi- 
pal component of the electric-held gradient tensor. The 
point quantity n zz is the main component of the tensor 
riij = (didj — y<5y)n(0), which has via the Laplacian 
two derivatives more than the main component of the 
EFG tensor V zz . n zz can be shown to be proportional to 
(y r 2m/i' 5 ) and therefore bears similarities with the elec- 
tric quadrupole field (Y 2m /r 3 ) as well as with the electric 
hexadecapole field (F 4m /r 5 ), cf. Eq. J7J. 



III. OBSERVABLE CONSEQUENCES 

All entries in the (classical) Tab. U correspond to an 
experimentally observable correction to the total energy. 
The first row lists energy corrections which are a product 
of scalar quantities. The leading term after the monopolc 
contribution MI (or -Eg") is the first order monopole shift 
MS*- 1 ', which experimentally manifests its presence in the 
well-known isomer. The second order monopole shift 
MS' 2 ) is only very rarely taken into account. One ex- 
ample where it matters is the case of muonic atoms 0, H[ 



7=2 
1 2 



mi = ±-, 



2v QI 



VQI 



2(v Q i + v QS ) 



VQI + VQS 



mi 



±1 



■} ■ 



Quadrupole + Quadrupole 
Interaction Shift 

FIG. 1: Energy levels for a nuclear spin of I = 5/2. This 
picture is not on scale: the shift of the levels as indicated by 
the arrow is in the most favorable cases (=heavy nuclei) only 
0.1 % of vqi. 



(atoms where a muon rather than an electron orbits the 
nucleus). Because a muon is much heavier than an elec- 
tron, its orbit is much smaller and the overlap with the 
nuclear charge distribution becomes much larger. This 
makes the second order monopole shift for muons much 
larger than it is for electrons [97| . 

All entries in the second row of Tab. [J are dot prod- 
ucts between spherical tensors of rank 2. The first one 
is the quadrupole interaction term QI, which splits ac- 
cording to Eq. (|2"2"1) energy levels that were degenerate 
under the monopole term. An example for the axially 
symmetric case (77 = 0) and nuclear spin / = 5/2 is 
given in Fig. [TJ The second term in the second row is 
the first order quadrupole shift QS^\ which shifts the 
energy levels that were split by the quadrupole Hamilto- 
nian, but preserves its overall symmetry (Fig.[TJ example 
for 7} = t]qs = [2H ) . The frequencies (energies) that set 
the scale for the quadrupole and quadrupole shift split- 
ting are (still considering 77 = r/gs = 0): 



VQS 



eQV zz 
h 

eQn zz 
Ue h 



(32) 
(33) 



(For the sake of lighter notation, we will use from here on 
vqs rather than Vq S (V) : we will not consider second order 
quadrupole shifts and therefore no confusion will be pos- 
sible.) The quadrupole shift does not change the overall 
symmetry, which in the example of Fig. Q] means that the 
1:2 ratio between the two energy differences is preserved. 
An experiment that measures such energy differences is 
not able to distinguish between the contribution by vqi 
and the one by vq$: it measures their sum only. A dis- 
cussion of the trends in the order of magnitude of the 
quadrupole shift will be given in Sec. IV CI and several 
experimental and computational strategies to exploit the 
quadrupole shift will be suggested in Sec. IVII 

Finally, the third row in Tab. [T] lists dot products be- 
tween tensors of rank 4. The leading term here is the 
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hexadecapole interaction for point nuclei. This term 
can in principle be distinguished experimentally from a 
quadrupole interaction because its symmetry is different 
(for instance, in Fig. Q] the 1:2 ratio would be slightly vio- 
lated). The HDI appears only for nuclei with I > 2, since 
only they have hexadecapole moments (21 > I rule for 
2 l multipole moments). Whereas the QI is well known 
and experimentally accessible by e.g. NMR or Molecular 
Beam Spectroscopy (see Sec. I VI A")) , the situation for the 
HDI is different. Since it was reported in 1955 for the 
first time [l2j , it has gone through cycles of confirmatory 
measurements and refutations. An overview is given in 
Ref. 0. 



IV. COMPUTATIONAL ASPECTS 
A. Formulation in spherical notation 

The electronic part n zz of the quadrupole shift will be 
calculated with a first-principles code and must therefore 
be translated in spherical form as it is common in such 
codes: 



0.5 



2 15 v 1 M 

T n ~~2 n 20(O- 



(34) 



The spherical component of the density, n2o (r) , which en- 
ters this expression, is the radial part of the (1 = 2, m = 0) 
component of expansion of the density n(r) in spherical 
harmonics: 



n (r) =y^n lm (r)Yi m (iY). 



(35) 



i in 



The I = 2 components are closely related to Cartesian sec- 
ond derivatives [l3T ] , which is the reason why they appear 
in the electric-field gradient and related quantities. 



B. Computational details: the FPLO code 

All calculations in this paper have been performed 
by the Density Functional Theory solid state code 
FPLO [II (version 8.00-31), which is is a full-potential 
band structure scheme and based on linear combina- 
tions of overlapping non-orthogonal atom-centred or- 
bitals. The core relaxation is properly taken into ac- 
count (so called all-electron method). We used the 
Local Density Approximation (LDA) for the exchange- 
correlation functional [l5j. FPLO can perform non- 
relativistic, scalar-relativistic as well as fully-relativistic 
calculations QjJ. In the latter, the Dirac Hamilto- 
nian with a general potential is solved. Recently, a finite 
nucleus has been implemented in FPLO, which is crucial 
for the present work (Sec. HVCj) . 
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FIG. 2: The density component ri2m(r) including inset, which 
zooms in the region around r = ) for a point nucleus (p.n., 
dashed lines) and a finite nucleus (f.n., full lines) plotted in de- 
pendence of r. The different methods are indicated by differ- 
ent colors: non-relativistic (NREL, yellow), scalar relativistic 
(SREL, red) and full relativistic (FREL, green). This calcu- 
lation for the hep metal Re was done by FPLO. All quantities 
are given in atomic units. 



C. Relativity and the role of a finite nucleus 

In order to obtain n zz , the limit of n2o (r) / r 2 for r — > 
must be calculated, cf. Eq. I|34p. It matters whether 
this is done within a non-relativistic (NREL), a scalar- 
relativistic (SREL) or a folly relativistic (FREL) frame- 
work. In the NREL or FREL formulations (no matter if a 
point or a finite nucleus is used in the calculation) ri2m (0) 
is exactly zero as it should be due to angular selection 
rules (Fig. [2]). In the SREL approximation, the (1 = 2, m) 
density, created from two divergent V\/2 functions, is to 
some extent wrongly non-zero at r = 0. This makes 
SREL-based methods (with or without a point nucleus) 
essentially useless for calculating properties that depend 
on n2a(r — >0), and we will therefore not consider SREL 
any further. 

For a point nucleus, the ratio of n2o(r) and r 2 con- 
verges for the limit r— >0 in a NREL formulation, but not 
in a FREL formulation (Fig. [3]). Since this ratio at r = 
is an observable quantity (see Eqs. (f3~4"]) and (|3"3"|) h the di- 
vergence for the better method (FREL vs. NREL) cannot 
be physical. And indeed, the divergence disappears if the 
approximation of a point nucleus is dropped and a finite 
nucleus is used in the calculation (Fig. [3]). Numerical 
values for this ratio turn out to be much larger for FREL 
compared to NREL, especially for heavy elements. 
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FIG. 3: The density component ri2 m (r)/r 2 for a point nucleus 
(p.n., dashed lines) and a finite nucleus (f.n., full lines) plot- 
ted in dependence of r. The different methods are indicated 
by different colors: non-relativistic (NREL, yellow), scalar 
relativistic (SREL, red) and full relativistic (FREL, green). 
This calculation for the hep metal Re was done by FPLO. All 
quantities are given in atomic units. 



The divergence of n zz in a fully relativistic point nu- 
cleus calculation might appear to be worrying at first 
sight. Wouldn't that mean that the quadrupole shift in 
Eq. (|33p is infinite? The answer is: no, because the oper- 
ator corresponding to n zz (Eq. (|3Tj) ) docs not have to be 
evaluated in the ground state for the point nucleus (which 
is the case that diverges at r — 0), but in the ground state 
after having added the two perturbations of Eq. (|20[) that 
describe the effect of quadrupolarly deformed finite nu- 
cleus (where the divergence is absent). The latter ground 
state can be constructed from the ground and excited 
states of the point nucleus case, applying the common ex- 
pression for the eigenfunctions in first order perturbation. 
This would, however, lead to rather lengthy expressions 
and to the inconvenience of having to use excited states. 
A pragmatic workaround is to use instead the ground 
state as calculated in a first-principles code that takes a 
finite nucleus into account. This is hardly an approxima- 
tion, as it was exactly the purpose of the perturbations 
in Eq. (|20[) to express the presence of a finite nucleus. 
Therefore, we conclude that the quadrupole shift can be 
obtained by evaluating the operator for n zz in Eq. $31]) 
for the ground state of the atom, molecule or solid calcu- 
lated fully relativistically and with a finite nucleus taken 
into account. This quadrupole shift has to be added to 
the contribution obtained by evaluating the operator for 
V zz in Eq. (|30|) in the ground state of the point nucleus 



case (and not in the ground state of the finite nucleus 
case, as the regular QI is really a perturbation to the 
point nucleus). 



D. Comparison with the PCNQM method 

We have described in the previous sections a proce- 
dure to obtain the influence on the quadrupole interac- 
tion of electron penetration in a finite nucleus by two 
subsequent applications of first order perturbation the- 
ory combined with finite nucleus calculations (Eq. ([55]) 
and Figs. and [3]). An alternative to this procedure 
is the point charge nuclear quadrupole moment method 
(PCNQM) plEal!]]: the electric-field gradient is not 
obtained as the expectation value of an operator, but is 
determined from the way how the total energy of the 
system changes upon inserting an artificial array of point 
charges around the nucleus. In this method, only total 
energies are required to obtain the clcctric-ficld gradient, 
which makes it particularly useful when the proper op- 
erator for the field gradient is not explicitly known. The 
latter is for instance the case as soon as a finite nucleus 
is used (Eq. (|2"2")) is valid for a point nucleus only), or 
for fully relativistic calculations at the 2-componcnt level 
(a complicated and not yet performed 'picture change' 
transformation would be needed to find the 2-component 
version of the EFG operator [HI.) The difference in 
EFGs between a 'finite nucleus + PCNQM' calculation 
and a point nucleus calculation (either with the regu- 
lar EFG operator or with PCNQM) gives the effect of 
electron penetration in the nucleus. One case where this 
difference is explicitly calculated is for 127 I in Lil (Ref. @ 
and Fig. SI). With PCNQM, the quadrupole shift can be 
obtained only numerically: there is no analytical expres- 
sion as Eq. (f3"3")) . 

In passing, we note here that a method that is quite 
analogous to PCNQM has been recently developed for 
the first order monopole shift correction MS'- 1 -' (isotope 
shift, isomer shift) [2l|, [22j as well. 



V. NUMBERS AND TRENDS 

In the present section, we will perform actual calcu- 
lations with the formalism described in Sees. [TTI and IIV1 
and examine trends in the relevant quantities: the nu- 
clear quasi-quadrupole moment Q, the electronic point 
property n zz , and their product: the quadrupole shift 
vqs- 



A. Trends in Q 

Consider a phenomenological model for a nucleus: a 
deformed sphere, with a radius R(0) given by [23j: 



R{9) =a(l + /3 2 Y 20 {9) + ^Y A0 (6) + ...). 



(36) 
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TABLE II: The nuclear radius a, quadrupole moment Q, de- 
formation parameter 02 and quasi quadrupole moment Q of 
a few isotopes. 



Isotope 


a [fm] 


Q [fm 2 ] 




Q [fm 4 ] 


9 Be 


2.84 


5.3 


0.22 


43 


47 T ; 


4.61 


30.2 


0.09 


644 


nl Cd 


5.95 


83.0 


0.07 


2 934 


138 La 


6.34 


45.0 


0.03 


1808 


179 ffi 


6.84 


379.3 


0.15 


17 760 


187 Re 


6.93 


207.0 


0.07 


9 945 


189 0s 


6.95 


85.6 


0.03 


4138 



where a is called the monopolc radius and the pi are 
deformation parameters. The monopolc radius depends 
in the first place on the atomic mass number A of the 
nucleus, and the main trend throu gh a lot of experimental 
values can be summarized by [24| |99j 

a(A) = 1.489 A - 294 fm. (37) 

Values for /3 2 fall rarely outside the range [—0.3, +0.3] 
(Ref. [25[ in combination with Eq. (|38p). As @4 is even 
smaller and enters only quadratically in the expressions 
we will need, it can be neglected for our purposes. Keep- 
ing only the linear order for /?2, we can now express the 
quadrupole moment and the quasi quadrupole moment 
in terms of a and P2' 

eQ ~ 3i/^/3 2 a 2 (38) 

V 5 Z7T 

eQ ~ a 2 • eQ. (39) 

The term quadratic in /3 2 as well as the quadratic $4 
term give corrections to Eqs. (|3"8"|) and (f3"9")) at the level 
of a few percent only, while they make the expressions 
considerably more involved - see Ref. @ . 

By Eqs. (|35|) - (|39"f . one can get a reasonable estimate 
for Q by inserting the monopole radius from Eq. (|37|) 
and the experimental quadrupole moment Q (e.g. from 
Ref. [H, Hil [27} ) ■ In this way, we obtain values for Q in 
the order of 10 4 — 10 5 fm 4 for heavy elements (Tab. HI)) . 
The Eqs. (|38[) - (|3"9"1) show that in order to get a large quasi 
quadrupole moment Q, the nucleus should be large (a is 
large) and strongly deformed (Q or /3 2 are large). The 
former implies heavy elements, while the latter is most 
easily fulfilled for heavy elements as well. 

B. Trends in n zz 

In order to get a feeling for the order of the magnitude 
of the electronic parts of the 0(2) interactions in Tab.fl] 




20 40 60 80 100 
-12- Z 
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A 

FIG. 4: The logarithm of the ratio of vqs and vqi as a func- 
tion of the mass number A. Blue diamonds: artificial crys- 
tal structures (see text and Ref. ||), fitted by the red line 
(Eq. (|40|) ). Orange triangles (down): experimental crystal 
structures (see Tab. IIII| ). The yellow triangle (up) [f|, green 
square |(| and red circles [l[ are values from the literature, 
see text. Inset: the same data but now as a function of Z, 
fit by Eq. l)4ip. The nuclear (green) and electronic (orange) 
contributions of Eq. I|42p are shown as well, shifted to match 
in the endpoint. 



we have calculated both V zz (the electronic part of the 
QI) and n zz (the electronic part of the first order QS) 
for some hexagonal close-packed (hep) metals throughout 
the periodic table. The results are shown in Tab. IIIII 
Both quantities increase with the mass of the element. 
But compared to V zz , which increases over two orders of 
magnitude, n zz is much more sensitive to the mass of the 
element and increases over eight orders of magnitude. 

In order to verify to which extent this conclusion ob- 
tained from Tab. IIIII is valid for other crystal struc- 
tures than hep, we investigated two series of purpose- 
built body-centered tetragonal (bet) crystals with differ- 
ent c/a ratios (0.8 and 1.2), and this for several elements 
throughout the periodic table. The results are reported 
in Ref. [9j and show the same trend as Tab. IIIII We con- 
clude that the mass of the element has a larger influence 
on the magnitude of n zz than the lattice parameters or 
the crystal structure. 



C. Trends in the quadrupole shift 

How do the nuclear part from Sec. IV Al and the 
electronic part from Sec. IV Bl combine to produce a 
quadrupole shift? The frequencies vqj (for the QI - 
Eq. ([32])) and v QS (for the QS - Eq. ^) for a set of 
hep and bet metals are reported in Tab. IIIII together 
with their ratio W os ly oi I ■ Experimental lattice parame- 
ters were used [28|, |29(, and n zz and V zz were determined 
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TABLE III: For a few atoms/nuclei that experimentally condense in the hep crystal structure (except for Pa, bet), this table 
lists the nuclear properties Q and Q (determined as in Table [IT}, the electronic properties V zz and n zz /eo (calculated by FPLO, 
see text), the quadrupole vqj and quadrupole shift vqs frequencies they give rise to (Eqs. (|32p and (J33J) (mind the different 
MHz and kHz units), and the ratio of the latter. 



Isotope 


1 


Q 


Q 


1 7 

Vzz 


n zz /eo 


VQI 


VQS 


\VQS/VQI\ 






[fm 2 ] 


llm J 


[1 n2 It' 1 2 1 

[10 V/m ] 


[10 V/m J 








9 Be 


3/2 


5 


42 


-0.08 


-6.07-10 -2 


-0.1 


10" 8 


5-10" 9 


47 T j 


5/2 


30 


644 


1.61 


3.27-10+ 3 


11.8 


-0.04 


3-10" 6 


49 T ; 


5/2 


25 


539 


1.61 


3.27-10+ 3 


9.6 


-0.04 


3-10" 6 


m Cd 


5/2 


83 


2 934 


7.48 


2.94-10+ 5 


150.0 


-14.9 


LHT 4 


177 Hf 


7/2 


337 


15 652 


7.89 


1.26-10+ 6 


642.3 


-341.4 


5-HT 4 


179 Hf 


9/2 


379 


17 760 


7.89 


1.26-10+ 6 


723.9 


-387.4 


5-10" 4 


185 Re 


5/2 


218 


10 386 


-5.51 


-1.8140+ 6 


-290.3 


324.9 


L10~ 3 


187 Re 


5/2 


207 


9 945 


-5.51 


-1.81-10+ 6 


-275.6 


311.1 


1-10" 3 


189 Os 


3/2 


86 


4138 


-6.65 


-2.9L10+ 6 


-137.6 


208.1 


2-10~ 3 


231 Pa 


3/2 


-172 


-9 357 


15.14 


8.1L10+ 6 


-629.8 


1309.7 


2-10" 3 



fully relativistically with a finite nucleus for n zz and a 
point nucleus for V zz (see Sec. HVBl) . Q was taken from 
the literature [27j and Q was determined as explained in 
Sec. IV Al The trends of n zz and Q to be larger for heavy 
elements, cooperate to produce a vq$ of which the rel- 
ative importance with respect to vqi is rather smoothly 
increasing with the atomic number A. 

This can be seen more clearly in Fig. |4] (blue dia- 
monds), which summarizes results for a larger set of 
28 elements in different crystal structures: hep with 
c/a=1.633 and 0.8 and bet with c/a=1.2 and 0.8, always 
with the experimental volume per atom (details are given 
in Ref. [§]). These data can be fit with the simple func- 
tions 

\v QS \ = 5.46 -Hr 12 ^ \v Q i\, (40) 
\v QS \ = 3.26 -10- U Z 4 \v QI \, (41) 

which are shown in Fig. 2] (red lines). The orange tri- 
angles (down) in Fig. 2] correspond to the experimental 
crystal structures from Tab. IIIII - they accurately follow 
the same trend. 

By taking the ratio of Eqs. (|3"3"|) and (|3"2")) and by 
filling out the lowest order expressions for Q and Q 
(Eq. p8]) ). the following simple analytic analogue for 
Eqs. (gnj) or (gU) is obtained: 

"- = (-5°^)"- < 42 » 

Since a = 1.26 Z 1 / 3 fm (obtained from the data of 
Ref. [13] plotted as a function of Z), the nuclear part 
a 2 scales with Z 2 / 3 . In order to fulfill the observed 
Z A dependence in Eq. (|4"Tj) . the electronic part should 



scale with Z 10 / 3 : n zz /(e V zz ) = 2.87 • 10~ 10 Z 10 / 3 far 2 . 
These two contributions are shown as the green (nu- 
clear) and orange (electronic) lines in the inset of Fig. 2] 
From this picture, it is clear that the electronic term 
contributes most to the increase of the quadrupole shift 
with A or Z. From Tab. IIIII we see that this is due to 
the strong increase of n zz . 

Eq. (|4"Tj) provides a quick way to estimate the order 
of magnitude of the quadrupole shift, for any element in 
any crystal structure, and without the need for a finite 
nucleus calculation. The only quantity that is required 
is vqi, which can be provided by several first-principles 
codes. As the scatter of the data points for heavier ele- 
ments shows, such an estimate can be one order of mag- 
nitude above or below the actual value. For isotopes with 
A > 175 (Z > 60), the quadrupole shift can reach 0.1- 
1.0% of the regular quadrupole interaction. 

There arc a few cases reported in the literature from 
which QS information can be deduced. These arc shown 
in Fig. 2] as well. The yellow triangle (up) was calculated 
by J. Thyssen et al. with a method very similar to ours 
for the single case of the Lil molecule. They found the 
ratio \vqs/vqi\ for 127 I to be 5 • 10~ 5 . The green square 
corresponds to 127 I in the same Lil molecule, obtained by 
the PCNQM method by Van Stralen and Visscher 0. A 
few estimates for the quadrupole shift obtained by sec- 
ond order perturbation theory were published in 1970 by 
P. Pyykko [l|. Those estimates were given relative to a 
pseudo quadrupole interaction only (|A II and Ref. [30l]). 
After converting these numbers, it turns out that for the 
LiBr molecule the ratio of vqs and vqj is about 10~ 10 
for 6 Li and 10 -6 for 81 Br (red circles in Fig. g]). These 
numbers follow the same trend as the quadrupole shift in 
first order perturbation, but are 1-2 orders of magnitude 
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smaller - this might be due to the fact that these were 
non-relativistic calculations. 



D. Other small perturbations to the quadrupole 
interaction 

When dealing with a quadrupolc-likc interaction that 
is as small as the quadrupole shift, it becomes relevant 
to take into account similarly small quadrupole- like inter- 
actions and perturbations of the quadrupole interaction 
that have a different origin. Some of these interactions 
have been proposed decades ago, in an era during which 
it was impossible to compute accurate values for them. 
Given the enormous advances in the possibilities of first- 
principles calculations since that time, it is worthwhile 
to list these effects here, to discuss them shortly, to put 
them into a general picture and to refer to the original lit- 
erature. This is done in App. El where we will deal with 
second order effects of magnetic origin, the isotopologuc 
anomaly and the influence of temperature. 



VI. EXPERIMENTAL IMPLICATIONS OF THE 
QUADRUPOLE SHIFT 

A. Accuracy of quadrupole interaction 
experiments and calculations 

In order to see whether or not the presence of the 
quadrupole shift can be experimentally detected, we 
should assess the accuracy that can be achieved in 
quadrupole interaction experiments. In order to find out 
which kind of information can be extracted if such ex- 
periments are combined with first-principles calculations, 
the best achievable accuracy in such calculations will be 
discussed as well. 

Experimental methods can be either non-radioactive 
ones as nuclear magnetic resonance (NMR) spectroscopy, 
nuclear quadrupole resonance (NQR) spectroscopy, laser 
spectroscopy (LS) and molecular beam spectroscopy 
(MBS), or radioactive ones as Mossbauer spectroscopy 
and perturbed angular correlation (PAC) spectroscopy. 
These methods can be applied to atoms and molecules 
(LS, MBS, NMR, NQR) or to solids (NMR, NQR, 
Mossbauer, PAC). A typical NQCC vq is of the order 
of magnitude of 100 MHz. These are the lowest achiev- 
able experimental error bars on vq for each method: 
5 kHz for NMR or NQR on single crystals with an axi- 
ally symmetric EFG [fEHl, 100 kHz for NMR or NQR 
on powder samples with a non-axially symmetric EFG 
[HI, [Hj|, 5 MHz for LS on atomic beams_|3|L 5-20 Hz 
(!) for MBS (H HI, 500 kHz for PAC lUl HI and 
Mossbauer spectroscopy [12 . When compared to the the 
quadrupole shift values in Tab. IIIII (typically 100 kHz 
for A>150), it is clear that only NMR in solids and es- 
pecially Molecular Beam Spectroscopy on molecules are 



sensitive to the quadrupole shift - provided the isotope 
under consideration is sufficiently heavy. 

First-principles calculations in solids are commonly 
done at the level of density functional theory (DFT), or 
with DFT as a starting point. DFT has been used with 
considerable success to calculate electric-fiel d g radients 
in solids, see e.g. Refs. [H, EE EH, E3, HI H EE EE 
E3, El, El, HE M, H2- As a rule of thumb, the DFT 
prediction is within 10% of the experimental value. 

First-principles calculations for (small) molecules can 
resort to Hartree-Fock calculations with correlation cor- 
rections. These are computationally much more demand- 
ing, but can in principle provide an arbitrary high preci- 
sion. The recent literature [H, [H, HE HI, E3, EH shows 
that correlation corrections using coupled cluster theory 
with single, double and (perturbatively) triple excitations 
(CCSD(T)), combined with sufficiently large basis sets 
and - where needed - with a (semi-)relativistic Hamilto- 
nian, provides highly accurate EFGs for small molecules. 
It has been claimed [12 that in this way an absolute ac- 
curacy with 4 significant digits can be reached. This is 
considerably better than the accuracy which DFT can 
provide for the EFG in solids. For molecules that are 
too demanding for a CCSD(T) treatment, DFT with the 
recently proposed CAMB3LYP* functional can be an al- 
ternative [6(|. DFT for EFGs in small molecules can be 
very unreliable [12 ■ 

At non-zero temperatures, vibrational states will be 
populated in solids and molecules, and in molecules ro- 
tational states as well. This will influence the electric- 
field gradient. In solid state calculations, this has so far 
only rarely been taken into account [olj |. In molecules, 
the effect of temperature is routinely taken into account 
in calculations [50, [6211 as well as in the analysis of ex- 
periments [H, l64l . |65| |. This allows an even more de- 
tailed comparison between experiment and theory for 
molecules. 



B. Determination of Q and Q: method 

With the experimental accuracies listed in the previous 
section, it is clear that experimental nuclear quadrupole 
coupling constants vq for NMR on single crystals and for 
molecular beam spectroscopy on molecules are affected 
by the quadrupole shift. This means that the experi- 
mentally determined value for vq would have a different 
value (outside the error bar) if the quadrupole shift could 
be "switched off". It does not mean, however, that by 
such an experiment the quadrupole shift itself can be de- 
termined: the QS manifests itself as an addition to the 
regular quadrupole interaction, and is indistinguishable 
from it (Eqs. (32J and ([55])): 

VQ « VQI + VQS = Q^-Q^- ( 4 3) 

The second term of this equation is even in the most fa- 
vorable cases 2-3 orders of magnitude smaller than the 
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first term (Fig. |4j). If vq is measured, if V zz is calculated 
from first-principles and if vqs is neglected, then the 
quadrupole moment Q can be determined from Eq. (j43]l . 
This has become the preferred procedure to determine 
nuclear quadrupole moments (e.g. [27l. l27l. I-431. loOl. lo6l. l67l. 

HEl). 

If V zz could be calculated with an arbitrary high pre- 
cision, the precision of the resulting Q is limited by ne- 
glecting vqs- One could choose not to neglect vqs, and 
apply Eq. (|43|) to at least two vq measurements in order 
to determine simultaneously a more precise value of Q 
and Q (or Q and a 2 ). This would be meaningful only in 
cases where the absolute deviations on the computed V zz 
and n zz values are small enough to make the uncertainty 
in vqi smaller than the value of vqs- The only hope to 
realize this is in the case of sufficiently heavy elements, 
for which, however, it might not yet be feasible to achieve 
the requested computational accuracy. 



C. Quadrupole moment ratios: the quadrupole 
anomaly 

When it is not possible to know experimentally the 
value of a quadrupole moment with sufficient accuracy, 
the next best thing to know are ratios of quadrupole mo- 
ments for two different isotopes, or for two different iso- 
meric states of the same isotope. As soon as a later 
experiment succeeds to determine one of the quadrupole 
moments in the ratio, the other one is known as well. 

The ratio Q1/Q2 of two quadrupole moments is com- 
monly measured as the ratio vq,i/vq,2 of two nuclear 
quadrupole coupling constants. Indeed, in the absence of 
a quadrupole shift, both ratios are identical if the two iso- 
topes or isomers are in the same environment and there- 
fore experience the same V zz (Eq. (|43[) h The presence 
of the quadrupole shift, however, spoils the equality of 
both ratios. It is straightforward to show that the ratio 
of quadrupole coupling constants is equal to 



with 8 



Hi 



14e y, ; 



O af 



(44) 



This formulation is strongly reminiscent to the Bohr- 
Weisskopf effect [7(| for magnetic hyperfine interactions, 
where the ratio between two magnetic hyperfine inter- 
action frequencies for two isotopes/isomers at identical 
sites is given by 



P-2 



Hi 

/'2 



(1 + A). 



(45) 



Here \i\ and [ii are the nuclear magnetic moments of the 
two isotopes/isomers, and A is the hyperfine anomaly. 
The ratio can be determined from hyperfine ex- 

periments on the two free isotopes/isomers in a known 
externally applied magnetic field. Comparison with 
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FIG. 5: The logarithm of the quadrupole anomaly 8 as a 
function of the mass number A, as given by Eq. (|46p . The 
value of n indicates the mass number difference between the 
heaviest (A) and lightest (A — n) isotope. The curves for 
n = 2 (full red line), n — 5 (blue dashed line) and n = 10 
(yellow dot-dashed line) are shown. For 47 Ti, 179 Hf and 187 Re 
(all n = 2) log \8\ can be calculated explicitly from Tab. Mil 
(orange triangles). 



the ratio as determined from experiments with the iso- 
topes/isomers incorporated in solids or molecules pro- 
vides the value for A, which can be as large as 2 % for 
heavy elements like 185 187 Re [tTJ] . A is nonzero because 
electrons that penetrate the nucleus do not interact with 
a point nucleus magnetic moment but with the spatial 
distribution of the magnetic moment over the nuclear 
volume. This slightly affects the effective hyperfine field. 
Therefore, the hyperfine anomaly is sensitive to the de- 
tails of nuclear structure, and can be used to test theo- 
retical nuclear models. 

In the same way the 8 from Eq. (|4"4")l - which can be 
called in analogy the quadrupole anomaly - probes de- 
tails of the nuclear charge distribution by electrons that 
penetrate into the nuclear volume. From Eq. (|4"4"|) , it 
can be seen that 8 is sensitive to the electronic quanti- 
ties, n zz and V zz , and the difference between the squared 
monopole radii of the two isotopes/isomers that are in- 
volved. 

In order to find a general trend and order of magni- 
tude estimate for 8, we combine the analytical function 
of Eq. d42]) with the fitted function of Eq. (gO]) and the 
square of Eq. (|37[) to obtain a numerical approximation 
for the electronic part n zz / {lAeoV zz ) in Eq. (|4"4")l . By in- 
serting this and the square of Eq. ([37)) for two different 
isotopes in the definition of 8, the following dependence 
of \8\ on the isotope mass number emerges: 



\S(A)\ 



5.46 



10 -12 A 3.079 



(A 



(A - n) 



)■ 

(46) 
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TABLE IV: Ratios of experimental quadrupole coupling constants for two different isotopes in two different diatomic molecules, 
collected from the literature. Only cases where the error bar on this ratio has been determined directly from the fit to the 
experimental data are reported (this error bar can be slightly different from what one would obtain using the error bars on the 
individual frequencies, see the discussion in Ref. [35|). The experimental value of the EFG (in 10 21 V/m 2 ) and the estimated 
value of 8 (Eq. (|46[) ~1 are given as well. 



Molecules Isotopes v na /v ma V^ v |<5| Ref. 

6 Li 19 F, 7 Li 19 F 6 Li/ 7 Li 0.020161 ± 0.000013 -0.44 5.9 • 10~ 10 [72] 

6 Li 127 I, 7 Li 127 I 6 Li/ 7 Li 0.02028 ± 0.00014 -0.18 [63] 

4i K i9 F> 39 K i9 F 41 K/ 39 K 1.217699 ± 0.000055 -5.6 1.3 • 10" 7 [73] 

41 K 127 I, 39 K 127 I 41 K/ 39 K 1.2174935 ± 0.0000099 -3.0 [74J 

87 Rb 19 F, 85 Rb 19 F 87 Rb/ 85 Rb 0.4838301 ± 0.0000018 -10.7 9.6 • 10" 7 [35] 

S7 Rb 35 Cl, 85 Rb 35 Cl 87 Rb/ 85 Rb 0.483837 ± 0.000022 -8.2 [75J 



This expression estimates the order of magnitude of 8 for 
two isotopes with mass numbers A and A — n. Curves 
for log |<5(A)| for n = 2, 5 and 10 are shown in Fig. 03 We 
observe that the quadrupole anomaly strongly increases 
with A (or Z), due to the increase of n zz . Mass number 
differences of 10 yield a value for 8 that is an order of 
magnitude larger than mass number differences of 2. For 
the 3 elements in Tab. IIIII for which information for 2 
isotopes is provided, Eq. (|46|) can be compared by values 
obtained by filling out the quantities of Tab. IIIII directly 
into Eq. ([4"4"|) . The values are shown by the orange trian- 
gles in Fig. [5] and correspond to the red fit (n = 2). This 
comparison shows that Eq. (|46[) is within one order of 
magnitude indeed a good estimate for S. The experimen- 
tally achievable accuracy of quadrupole moment ratios is 
of the order of 10~ 6 (see Tab. HVj). This means that for 
many isotopes the presence of 8 affects the experimental 
values. 

Unfortunately, whereas in Bohr-Weisskopf experiments 
the unperturbed ratio can be determined from ex- 

periments on free nuclei in an externally applied mag- 
netic field, this is not possible for quadrupole interaction 
measurements: electric-field gradients that can be gener- 
ated by man-made devices are too small to allow mean- 
ingful quadrupole interaction measurements [76| . There- 
fore, a slightly different method has to be used. One 
could perform 4 quadrupole interaction experiments on 
two isotopes ('m' and 'n') of the same element, each of 
them being part of two different molecules ('A' and '£?'). 
For instance, m X in m XA and m XB molecules, and n X 
in n XA and n XB molecules. This yields four exper- 
imental frequencies z^, and u% . By applying 
Eq. (|4"4")l twice, it can be seen that the NQCC ratios are 
not necessarily identical to each other for the two differ- 
ent molecules, with the difference being determined by 

11>zz/Vzz- 




v m Qm ( 1 nf z 1 2 2 \ \ (a q\ 

^ = o:{ 1+ T^s^- a ^)- (48) 

As long as the quadrupole shift (cx n zz ) does not play 
a significant role, the two experimental frequency ratios 
at the left-hand side are within their error bars identical 
to each other. If, however, the quadrupole shift would 
be large enough, these two experimental frequency ratios 
would differ from each other. This is a completely experi- 
mental procedure to detect the presence of the quadrupole 
shift effect. Tab. IIVI lists a collection of experimental 
NQCC-ratios in diatomic molecules determined for three 
such sets of 4 experiments, which gives an impression of 
the experimental accuracy that can be achieved. The es- 
timated order of magnitude for |<5| (Eq. (|4*fl)) ) is given too. 
For none of these cases, 8 is expected to be large enough 
to affect the experimental ratios. Tab. IIVI combined with 
Fig. [S] suggests that if the best experimental accuracies 
of 10~ 6 can be achieved for isotopes with A>150, then 
the influence of 8 could be observed. The heavier the ele- 
ment and the larger the size-differences between the two 
isotopes, the more likely large (5-values are. Interestingly 
enough, the quadrupole coupling constant ratios for the 
two K isotopes in the KF and KI molecules differ from 
each other in the 4 th digit, and this difference is an or- 
der of magnitude larger than the experimental error bars. 
Given the estimate for 8, the quadrupole shift is expected 
to give an effect in the 7 th digit at best. It is therefore 
unlikely that this set of K-expcriments represents an ex- 
perimental observation of the quadrupole shift (it could 
be due to one of the other effects discussed in App. [5] or 
due to an experimental problem). Nevertheless, it would 
be interesting to perform similar experiments with the 
same accuracy for heavier elements, where 8 is expected 
to be larger. 

One step further is to solve the system of the two equa- 
tions (|4T|) and (|4"5)) for the unknown quantities Q m /Qn 
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and (a? 



v nb 14e V 2 a s 




n zz 




14e V z % 


14e V^ 




_ "ml 




v n b 






Vnb 14e V^ 





(49) 
(50) 



All quantities at the right-hand side of these equations 
can either be measured or calculated, such that the quan- 
tities at the left-hand side are effectively determined by a 
combination of experiment and theory. Clearly, this is a 
game with very small numbers. The difference between 
the two frequency ratios in the numerator of Eq. (|50|) 
is of the same order of magnitude as the S in Eq. (|44|) : 
10 -5 for heavy elements. The same considerations as 
in Sec. IVf Bl apply here: an extreme accuracy in experi- 
ments as well as in calculations is needed in order to get 
to a reliable conclusion. Furthermore, the procedure as 
described here can be disturbed by the presence of a few 
other small quadrupolc-likc effects that are discussed in 
Sec. |VD] and App. [ft] 



VII. CONCLUSIONS AND OUTLOOK 

In this work, we described how electron penetration 
in the nuclear volume leads to the quadrupole shift: a 
small perturbation of the regular quadrupole interaction, 
which depends on the second derivative of the electron 
charge density at the nucleus (n zz ), as well as on the 
size and shape of the nucleus (Q). An explicit expres- 
sion for the quadrupole shift that can be implemented 
in a band structure code was derived, and DFT calcu- 
lations were performed for a set of crystalline materials. 
It was shown that meaningful numerical values for the 
quadrupole shift can be obtained only for fully rclativis- 
tic calculations that take a finite nucleus into account. 
Therefore, the quadrupole shift is one of the few cases 
where the commonly used scalar-relativistic approxima- 
tion is definitely insufficient. 

The quadrupole shift is a small effect. Its order of 
magnitude appears to be related in the first place to the 
atomic number A of the clement under consideration, and 
to a lesser extent to the crystal structure (Fig. [4} . This 
is predominantly due to the way how n zz depends on Z. 
The quadrupole shift is orders of magnitude smaller than 
the regular quadrupole interaction for most elements, but 
can reach 1% to perhaps 10% near the actinide region. 

We have pointed out how the quadrupole shift can play 
a role in a more accurate determination of quadrupole 
moments and quadrupole moment ratios. The compari- 
son of two accurately measured quadrupole coupling con- 
stant ratios provides a purely experimental way to ob- 
serve the presence of the quadrupole shift. For suit- 
able cases, the required experimental accuracy can be 
reached by e.g. molecular beam spectroscopy. With fur- 



ther advances in the absolute accuracy of ab initio calcu- 
lations for n zz and V zz , awareness of the existence of the 
quadrupole shift will help to extract more precise nuclear 
information from quadrupole coupling experiments. 

Suggestions for further work are at the conceptual, 
computational as well as on the experimental level. Con- 
ceptual: it remains to be understood which features of 
the electron density are responsible for the observed in- 
dependence of n zz and for the dependence of n zz for a 
given element on the crystal structure. Understanding 
those mechanisms would help to single out situations 
where the quadrupole shift is maximized. Computa- 
tional: in the present work, only DFT calculations for 
solids were performed, whereas the most accurate exper- 
iments are available for molecules. DFT for molecules is 
not likely to provide very accurate results, but quantum 
chemical calculations can do much better in this respect. 
It would be interesting to examine for instance the value 
of the quadrupole shift for heavy elements in a set of 
molecules. Experimental: sets of 4 quadrupole coupling 
experiments as in Tab. IIV1 done for heavy elements and 
with high accuracy, provide a way to observe the pres- 
ence of the quadrupole shift experimentally. It would 
be most efficient to make first a computational study, to 
identify among those molecules that are experimentally 
most easily accessible the ones where a large quadrupole 
shift is most likely. 
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APPENDIX A: OTHER SMALL 
PERTURBATIONS TO THE QUADRUPOLE 
INTERACTION 



As announced in Scc. lVDl this appendix discusses sev- 
eral other small quadrupole-like effects that might be of 
similar magnitude of the quadrupole shift. 



1. Second order effects of magnetic origin 

Van Vleck, Rabi, Foley and Ramsey [z3, [zl, [79[ dis- 
cussed half a century ago a pseudo quadrupole interac- 
tion in molecules that has a magnetic origin. This has 
been later elaborated upon by P. Pyykko, especially for 
the case of metals [3(| H3] • Mathematically, this pseudo 
quadrupole interaction arises in the same way as the 
quadrupole shift when the latter is derived by second 
order perturbation theory (see the end of Sec. Ill C[) . 
The small perturbing Hamiltonians are now the ones 
that give rise to the magnetic hyperfine field: the nu- 
clear spin/electron orbit Hamiltonian (Hi), the nuclear 
spin/electron spin dipolc-dipolc Hamiltonian (H.2) and 
the Fermi contact Hamiltonian (Hz). In second order 
perturbation, the square of Hi, the square of H2 and 
the cross-term of H2 and H3 contain an I 2 operator that 
gives rise to a quadrupole-like interaction (compare to 
Eq. (|22[) . provided axial symmetry (77 = 0) is present). 
Such a term is included even in first order in Hi. These 
pseudo quadrupole interactions were shown to be at the 
level of 10 -4 -10 -6 of the regular quadrupole interaction 
in molecules (3(| 
to 1% in metals 



781 [79J , and reach in favorable cases up 
30[. The values in molecules are there- 



fore of the same order of magnitude as the quadrupole 
shift (Fig. Hand Tab. EI]). 

Strictly spoken, these quadrupole-like contributions 
have a different status than the quadrupole shift in 
Sec. |TTJ The quadrupole shift Hamiltonian (Eq. ([231) ) 
has exactly the same structure as the quadrupole in- 
teraction Hamiltonian (Eq. (|22|) ). and they are there- 
fore completely indistinguishable from each other. The 
quadrupole-like interactions discussed in the present sec- 
tion have in the first place a I 2 dependence which splits 
the nuclear levels in a quadrupole-like manner as long as 
the environment has axial symmetry (77 = 0). The im- 
portant case of linear molecules has this symmetry. In 
less symmetrical environments, one could in principle dis- 
tinguish between these quadrupole-like interactions and 
the quadrupole shift. As this symmetry breaking is it- 
self a small effect, however, such considerations are not 
expected to be of much practical value. 

For completeness, we mention two other sources of 
quadrupole-like interactions, which are believed [8(| to 
be even smaller than the previously described ones: the 
influence of an external magnetic field [HI, HH and nu- 
clear polarization [&|. 



2. The isotopologue anomaly 

Through high-precision molecular beam experiments, 
Cederberg et al. have drawn attention to the fact that 
the quadrupole interaction at nucleus B in an AB di- 
atomic molecule slightly depends on whic h iso tope is 
taken for element A. For the isotopologues |l00l | 7 Li 127 I 
and 6 Li 127 I, this isotopologue anomaly is 0.007% of the 
regular quadrupole interaction at 127 I (an absolute shift 
of 14 kHz) [H. For 41 K 127 I and 39 K 127 I, the relative ef- 
fect at 127 I was 10 times smaller (0.0007%, absolute shift 
of 0.6 kHz) [l3, while for 39 K 81 Br and 39 K 79 Br there 
was no effect found at all on 39 K [lij]. The orig in of 
the isotopologue anomaly is not understood [63|, but 
from the literature review we present in Tab. [Vj a cor- 
relation between the relative value of the isotopologue 
anomaly and the relative mass number change for the A- 
isotopc is clearly present. The 50% relative mass change 
between hydrogen and deuterium results in an isotopo- 
logue anomaly of 0.2-1.0%. It tends to be lower for the 
33% mass change between deuterium and tritium (0.1- 
0.6%), although the error bars prevent unambiguous con- 
clusions. Much smaller but definitely non-zero frequency 
differences are observed for Lil and KI as well (0.007% 
and 0.001%). Tab. fVl suggests that these isotopologue 
anomalies tend to become undetectably small for rela- 
tive mass number changes below 5%. 

The isotopologue anomalies as presented in Tab. fVl 
were obtained as the difference between (^=0, J=0)- 
terms in the vibrational/rotational expansion of the 
quadrupole coupling (sec e.g. Eq. (19) in Ref. [HI). 
This lowest order term is not exactly equal to the static 
quadrupole coupling constant at the equilibrium inter- 
nuclear separation, due to the presence of an additional 
constant (the aB 2 term in Eq. (18) of Ref. [34[) which is 
mass-dependent and therefore isotope-dependent. This 
aB 2 term could therefore be an obvious candidate to ex- 
plain the observed frequency difference. However, Ceder- 
berg et al. have shown for CsF that aB 2 is negligi- 
ble [93| , because it is an order of magnitude smaller 
than the (v=2, J=0)-term, which itself is experimen- 
tally known to be small. We verified that the same 
argument holds true for Lil (Am re j = 14%) [34| and 
HI (Am r ei = 50%) [H, [9l|. Therefore, it is safe to con- 
clude that the larger as well as the smaller frequency 
differences in Tab. IVl are not significantly influenced by 
the aB 2 term and represent a real difference between 
two static quadrupole coupling constants, a difference of 
which the origin remains to be understood. 

Isotopologue anomalies in the kHz region can be of the 
same order of magnitude as the quadrupole shift vqs- 
Their existence puts further limitations on the numeri- 
cal information that can be extracted from a comparison 
of experimental quadrupole coupling constants and first- 
principles calculations. Indeed, the first-principles val- 
ues for V zz and n zz that appear for instance in Eqs. (|40|) 
and (|50p can only be calculated for specified elements in 
the molecule, not for the isotopes. On the other hand, 
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molecule A 


molecule B 


V Q A 


V Q,B 


Av 




Am re ( 


Ref. 






(MHz) 


(MHz) 


(kHz) 


(%) 


(%) 




'H 81 Br 


2 H 81 Br 


447.9(14) 


443.363(105) 


4537(1500)* 


1.023% 


50.0% 


[841 


1 H 79 Br 


2 H 79 Br 


535.4(14) 


530.648(74) 


4752(1500)* 


0.896% 


50.0% 


[841 


1 H 37 C1 


2 H 37 C1 


-53.436(95) 


-53.037(113) 


-399(200)* 


0.752% 


50.0% 


[841 


1 H 35 C1 


2 H 35 C1 


-67.800(95) 


-67.417(98) 


-383(200)* 


0.568% 


50.0% 


[841 


l H 127j 


2 H 127j 


-1828.059(51) 


-1823.226(54) 


-4833(100)* 


0.265% 


50.0% 


[85 861 


2 H 35 C1 


3 H 35 C1 


-67.417(98) 


-67.0(6) 


-417(700) 


0.622% 


33.3% 


[841 


2 H 79 Br 


3 H 79 Br 


530.648(74) 


530(2) 


648(2100) 


0.122% 


33.3% 


[841 


2 H 81 Br 


3 H 81 Br 


443.363(105) 


443(2) 


363(2100) 


0.082% 


33.3% 


[841 


2 H 37 C1 


3 H 37 C1 


-53.037(113) 


-53.0(6) 


-37(700) 


0.070% 


33.3% 


[841 


6 Li 127 I 


7 Li 127 I 


-194.33834(20) 


194.35241(20) 


14.07(40)* 


0.007% 


14.3% 


[631 


35 Cl 45 Sc 


37 Cl 45 Sc 


68.2067(29) 


68.2062(29) 


0.5(6.0) 


0.000% 


5.4% 


[871 


39j^127j 


41 K 127j 


-85.471138(7) 


-85.471721(12) 


0.583(20)* 


0.001% 


4.8% 


[74] 


63 CuOC- 127 I 


65 KOC- 127 I 


-593.465(9) 


-593.485(10) 


20(20) 


0.003% 


3.1% 


[88] 


79 Br 39 R 


81 Br 39 K 


-5.032957(9) 


-5.032957(9) 


0.000(20) 


0.000% 


2.5% 


[74] 


79 Br 45 Sc 


81 Br 45 Sc 


65.2558(32) 


65.2597(38) 


-3.9(7.0) 


0.000% 


2.5% 


[89] 



TABLE V: The isotopologue anomaly for a set of diatomic molecules. First two columns: the molecules with their isotopes 
- the isotope for which the NQCC is measured is put in bold. Second two columns: the NQCC in MHz. A^: the difference 
between the preceding two columns (kHz) - cases where the error bar allows to conclude the difference is not zero, are labeled 
by a "*"'. Av re i: relative frequency difference (%). Am re ;: relative change in atomic mass number for the neighboring isotope 
(%). 



the purely experimental determination of the presence 
of a quadrupole shift by 4 NQCC measurements as in 
Eqs. (|47|) and (|48|) is not disturbed by the isotopologue 
anomaly, as long as one makes sure that the isotopes for 
A and B in the ( m ' n )XA and ( m '")XB molecules remain 
identical in all 4 cases. 



3. Temperature and vibrations 

The entire discussion so far implicitly assumed static 
molecules or crystals (0 K and no zero point vibrations). 
At non-zero temperatures, vibrational states will be 



populated, and in molecules rotational states as well. 
These will influence the electric-field gradient and 
therefore the quadrupole interaction. The effect is in the 
range of 1-10%, and should therefore certainly be taken 
into account in high-precision studies. In molecules, 
this effect can be described with high accuracy using 
a Schlier-Dunham treatment [HI, [63|, HH, [65j . and 
quadrupole coupling experiments are routinely analyzed 
according to this formalism [H, [63|, [HI, [83, [H, [9ll . 



Similar studies in solids are rare - see e.g. Ref. [6l| for 
hcp-Cd, where a contribution of 1.6% due to zero-point 
vibrations was found. 
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